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Abstract 

The stress analysis of pressurized circumferential pipe weldments under 
steady state creep is considered. The creep response of the material is gov- 
erned by Norton's law. Numerical and analytical solutions are obtained by 
means of perturbation method, the unperturbed solution corresponds to the 
stress field in a homogeneous pipe. The correction terms are treated as stresses 
defined with the help of an auxiliary linear elastic problem. Exact expressions 
for jumps of hoop and radial stresses at the interface are obtained. The pro- 
posed technique essentially simplifies parametric analysis of multi-material 
components. 
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1 Introduction 

The welded pipelines subjected to high pressure and temperature are widely used 
in different branches of industry. Under such conditions, the creep and damage 
effects should be ta ken into account for accurate assurance of long-term reliability 
( Roche et al. . Il992h . Ap plication of computatio n al continuum cree p damage me- 



I Ito cne et all IIMMJ ). Applica tion ol computatio n al continuum cree p damage me- 
chanics (see, for example lAltenhadi et all l|20mh . IrTavhurstI fcOOIah ) 

coupled with 



increasing power of computers can accomplish this task. In recent years the finite 
element method has become the w idely accepted tool for the structural analysis in 



the creep range ([Havhurstl . l2001al ). A user defined creep material subroutine with 
appropriate constitutive and evolution equations can be developed and incorporated 
into the commercial finite element code to perform a numerical time step solution 
of creep and long term strength problems. 

On the other hand, in addition to more and more sophisticated numerical anal- 
ysis, simplified models of creep response are required. These models should provide 
a better intuitive insight into the problem and give a quantitative description of the 
solution. 
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The assessment of reliability of user-defined creep material subroutines and the 
choice of suitable numerical parameters like the element type, the mesh density, 
and time step control are complicated problems, particulary if studying creep of 
multi-material structures. Therefore, it is important to have reference solutions 
of benchmark problems. Such solutions should be obtained by use of alternative 
analytical or semi-analytical methods which do not require the spatial discretization 
techniques and allow for studying the behavior of stress and deformation gradients. 
The objective of this paper is to develop an alternative semi-analytical solutions 
to creep problems for multi-material pipe structures. Particularly we address the 
analysis of stress gradients in the local zones of material connections. 

To obtain a semi-analytical solution we shall make the following simplifications. 
We assume the idealized material behavior having the secondary creep stage only. In 
this case the steady state solution of creep in the pipe exists, for which the stresses do 
not depend on time. We assume that the difference between the material properties 
of constituents is not great. Particularly the difference between the minimum creep 
rates for the same stress level should not exceed the value of 2. 

The lifetime of a welded pipe under creep conditions is less than that of ho- 
mogeneous one. The effect of reliability reduction is of big interest, therefore 
large numbers of model problems were proposed. The most commonly used ap- 
proach sim ulates weldment as a region with non-un if ormly distributed materia l 
properties dBrowne et al. I (ll98ll).IColeman et all (Il985h. lHall and HavhurstJ (|l99lh . 



Perrin and Havhurstl (|l999j ). iHavhurst et all (|2001bl ). iHvde et all (|2003l )). 

Within the framework of this approach it is often necessary to consider a number 
of parameter distribution cases. The amount of problems to be solved increases with 
the number of changing material parameters and parametric analysis becomes very 
complicated. 

Drawing an analogy with some s i mple systems, for which an analy tical solution 



Drawing an analogy witn some s i mple systems, tor wnicn an analy tical solution 
is available fH vde et al. l (|l99rtl200nUNaumenko and Altenbachl ifeOOSn ) is useful for 
understanding how the parameter change can affect the solution. However this is 
not enough for a proper estimation of stress distribution. 

Two main types of constitutive equations are often used in weld modelling, 
namely, Norton's steady -state creep law and continuum damage mechanics equa- 
tions for tertiary creep ( Kachanov . 1986). Below only the constitutive equations 
of Norton's law are considered. Even by this assumption the structural response 
may be captured very well. Some special techniques are used to predict t he fai l- 
ure life more precis e ly us in g steady-sta t e solution 
Nikitenko and Zaev 



and Zaevl Jl979h. iHyde et all (Il99sl Il999h . IPerrin et all l|200m ). As 
Browne et all (|l98ll ). IColeman et al.1 (|l985l) for the particular pipe welds 



Leckie and Hayhurstl (|l974| ) 



shown in 

investigated, the steady-state analysis underestimates the failure time by about 20- 
40 percent, but predicts failure position quite well. 

To simplify the parametric analysis, we study a family of boundary-value prob- 
lems for mult i- component pipe creep depending on a small parameter s. When s = 0, 
the problem is reduced to the case of homogeneous pipe creep. T he correspondin g 
solution crfO) f or the steady-state stress distribution is well known ( Odqvistl ( 19741 ). 
Malinin ( 1981 ). Hvde et al. ( 19961 )) and considered to be the basic solution. In order 
to get the common solution cr(s) from this basic one, correction terms should be 
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Figure 1: System configuration (a) and boundary conditions (b) 



added. The equations for the correction terms are formally obtained as perturba- 
tions of a boundary- value problem with respect to s. These equations formulate 
a problem of linear elasticity with respect to linear-elastic solid with anisotropic 
elastic properties. 

The utility of used technique is guaranteed especially because the theory of 
linear elasticity is in a very satisfactory state of completion; every complicated case 
of parameter distribution can be treated in a routine manner as a combination of 
simple ones. 

Correction terms are obtained numerically for some problems with the help of 
the Ritz method. In some cases the simplicity of geometry enables us to construct 
an approximate analytical solution. An exact expression for stress jumps at the 
interface is obtained. Numerical solutions of the nonlinear problem are obtained 
with the help of ANSYS finite element code for comparison. 

Throughout light-face letters we denote scalars, the bald-face letters stand for 
tensors. The notation ( ) is used in vector-matrix form of constitutive equations to 
designate vectors cr and e of stress and strain components. 

2 Problem formulation: two-material model of 
narrow gap we Id me nt 

In this section we consider a two-material model only. It will be shown later that 
the solution for some multi-material models can be reduced to this case. 

The configuration analyzed is shown in Fig. ^ Assume that (r,z, 0) G O. x 
[0,27t), (r, z, 0) G D~x [0,2n), (r,z, 0) G D + x [0,2n) describe the volume occupied 
by the solid, by the weld metal, and by the parent material, respectively. Here 
Q = [ Ti , tJ x [0, H), Q + = [n, r D ] x [H, H], D. = Or U Q+. The basic equations of 
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the problem are given below. 
1. Equations of equilibrium 



V-cr = 0, (1) 

where <x is the stress tensor. In the volumetric forces are ignored. 

2. Strain-displacement relations 

e = 1 ( Vu + Vu T ) , (2) 

where £ is the linearized strain tensor and u is the displacement vector. 

3. T he govern i ng eq uations can be summarized as Norton's creep law (see, for 



example. lOdqvistl (197 



£ = As(ff vM ) n \ s = a- ^tr(<r)I, ct v m = : s, (3) 

where ( ) is the time derivative, s is the stress deviator, ct v m is the von Mises 
equivalent stress, I is the second rank unit tensor, and A, n are material constants. 
Note that Norton's creep law is often written as £ = |a s (ovm) 11-1 - It means 
A = |a in our notation. 

If we take into account that the problem is axisymmetric, we have in cylindrical 
coordinates the following equilibrium equations 

3o> 1 , 3o rz 3cr rz 1 , 3ff z 

or r oz or r oz 

strain-displacement relations 

3u r u r 3u z 1 / 3u r 3u z \ 

3r r 3z 2 \ 3z 3r / 

and boundary conditions 

<j r = — p, cr rz = at r = Ti; cr T = o rz = at r = r Q ; (6) 

cr TZ = 0, u z = at z = 0, z = H. 

Finally we define the distribution of the parameter A(r, z) in Norton's creep law 
Q as a piecewise constant function in £1 

A = A+ in Q+ A = A - in Or. (7) 

The stress field a, defined by (JTJ) — (JJJ), does not change if instead of A(r,z) 
we use AA(r, z), where A is any nonzero constant. Consequently, without loss of 
generality it can be assumed that 1 

A = A+ = 1 in + , A = A - = 1 — s in O . (8) 



1 Here s is not necessary a small parameter. In practice, A + and A might differ essentially. 
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If we eliminate displacements from the strain-displacement relations (0), we ob- 
tain compatibility equations 

d(e) = C 2 (e)=0, (9) 



Ci(e) = r(e T - 9 ^ 9 ), C 2 (e) 
or 



3 2 £ 3e z -3e TZ 

r — — H 2 . 

3z 2 3r 3z 



(10) 



We consider the weak form of compatib ility equations e xpressed by the equation 
of complementary virtual power principle ( Washizu . 1982) 



L(cr, s) < 6cr >= 0, VScr, 



L(cr, s) < 8a >= 



e(cr, s) : ScrdO. 



Here we use the brackets < • > to enclose the argument of a linear operator; e(<7, s) 
is the strain rate defined by (jHJ) and (JHJ); dO = rdrdz; 5cr is a virtual stress field 
that satisfy the equations of equilibrium (@J) and homogeneous boundary conditions 

cr r = cr TZ = at r = T|, r = r D ; cr rz = at z = 0, z = H. (12) 
In what follows we search function cr(s), such that 

L(ct(s),s) < 5cr >= 0, VScr. (13) 

3 Perturbation method 



3.1 Unperturbed solution: creep response of a homoge- 
neous pipe 



Now suppose that A + = A , i.e. s = 0. The problem (jTTj) is r educed to a one - 
dimensional, a n d the solution cr° of this problem is well known (jOdqvistl ()1974l ). 
Malininl (fjflisi. iHvde et al.1 ifjoOfih ) 



a T ° = a + a T r- 2/n , a° e = a + a e r- 2/n , cr° = a + a z r- 2/n , c° rz = 0, (14) 



2/n 

r 2/n 2/n ' 1 



2/n 2/n 
2/n 2/n 

r -r/ 



n-2 



n 



a T , a z 



TV 1 



n 



3.2 First perturbation 

Let s be a small parameter. Assume that 2 



tr s 



cr 



scr + o s 



(15) 



2 The .i ustification of p erturbation method could be performed by means of Implicit Function 
Theorem ()AntmarJ . ll995|) . Precise conditions justifying perturbation method should be formulated 
in an appropriate function space. 
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Here cr 1 = d °j; | s= o is the unknown derivative which must satisfy the equations of 
equilibrium and the homogeneous boundary conditions (fT2"j); o(s) is the little-o 
Landau symbol. From (JHJ) it follows that 

L(cr,s) < Scr >= (L°(<r) + sL 1 (cr)) < Scr > (16) 

with L°(cr) = L(<r,0), V[a) = Substituting and © in 0, we get 

dL° 



L°(cr )+s. 



CT o< cr 1 > +sL 1 (cr°) + o(s) ) < Scr >= 0. 



dcr 



Since L (cr ) < Scr >= and < Scr >—> as s — > 0, we have an equation for cr 



dL° 
dcr 



r o< cr 1 >< Scr >= -iJfcr ) < Scr >, VScr. 



(17) 



We will analyze more closely the linear operator \ a=a o< ■ > in the next section. 



4 Auxiliary problem 



4.1 Linear elastic material 



In the previous section it was shown that the correction term cr 1 can be found from 
the linear equation (|17|) . It is clear that 



dL° 
dcr 



i< cr 1 >< Scr >= 



^(ff°,0) < cx 1 >: ScrdO, 

Off 



where 



£(cr,0)=s(cr) (a vM (o")) n - 1 . 
Let us introduce a vector notation ( ) as 

cr = (cr r , ere, cr z , cjYz) T , e = (e T , £e, £z,2£ rz ) T . 
Substituting (|2T)J) in (fT£J) and differentiating ([THjl we obtain 
dL° 



da 



I<x=<tO< f 1 >< Scr >= 



cr 1 ) 7 C Scr dO. 



a 



Here we have introduced the compliance matrix as follows 



c=|i (< f°) 

Off 



n-2 



2/3+(n-1)/2 -1/3-(n-1)/2 
2/3 + (n-1)/2 

sijm. 



(19) 



(20) 



(21) 



(22) 



-1/3 
-1/3 
-1/3 
2 
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H 



Figure 2: Displacement jump 

Thus, the left-hand side of (I17j) can be treated as an internal complementary 
virtual work with respect to linear elastic solid with constitutive law 



(e T) £9, £ z , 2e TZ ) T — C(cr T , a Q , cr z , cr rz ) T . 
The eigenvalues of the compliance operator C are 

(A 1) A 2) A 3 ,A 4 ) = (v/3— r- 2/n )- 2 (0,2,2,n) 



TV 



(23) 



(24) 



The problem of steady-state creep is reduced to the elasticity problem for or- 
thotropic, incompressible, and inhomogeneous solid (J23|) . 



4.2 Displacement jump 

We now seek to convert the right-hand side of (|17|) into a surface integral through 
Gauss theorem. It can be proved that 



V(v°) < 8a >-- 



i{a°,0) : So- dCl 



3 ^ a T |a T 



n-1 



1 



rn' 



-5a rz rdr. (25) 



Q- 



{z=h} 



Let us show that (J25|) prescribes a jump of displacements at the interface {z = h) 
(see fig. EJ). Consider the principle of virtual complementary work for linear elastic 
solids CI and One gets 



o--) T C ba dO 



(u T 6cr TZ + u z 5ff z )rdr, 



(26) 



n- 



{z=h} 



6cr T = 6cr rz = at r = r-t, r = r Q ; 5a rz = at z = 0. 



(27) 
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'cr+r C Sct dO = 



(u+Sa TZ + u z Scrjrdr, 



{z=H} 



Suppose 



5cr r = 5cr rz = at r = r-t, r = r Q ; Scr TZ = at z = H. 



°£ - °"rz = °. a+ - cr = at z = h; 



u z — u z =0, u r - u r = 3 2 a T |a 



at z = h; 



then it can be shown that the solution cr 1 of (|17|) has the form 

cr + in D + 



O" 



a in O 



4.3 Stress jump 

We define a notation for jumps of field variables at the interface {z 



v = v — v 



Substituting (|3T|) for u T in (j5J), we get 



[e T ] = -c/r 2 , [e ] = c/r 2 , c = 3^2 1 a r |a r | n 



If we combine this with (J22J) . ()23|). and (J3Dj), we obtain 



v^arlarl _ 4/n 



n _ v^arlarl _ 



4/tl 



4.4 Case of multi-material structure 

In this subsection we consider m-material structure (m > 2) 

Q = Q,UQ 2 U...U a m , 

= [Ti,Tj X [Zj,Z j+1 ], zi = 0, z m+1 = H, 
A = Aj in Cl h ) = 1,...,m; A m = 1 , A m _, ^ A T 
There exist unique s and <Xj , j = 1 , . . . , m — 1 such that (fig. 0J) 

m-1 

A(z) = 1 - s Y_ a i H ( z i+i ~ z ). <*m-i = 1 • 
5=1 



Here Hfz) is the Heaviside function. 
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At = 1 -s(l +ai) A 2 = l-s A 3 = 1 z 
• • • • ► 

Z-i = Z2 z 3 z 4 = H 

Figure 3: Distribution of A in m-material structure (m = 3) 

If we use the small parameter method to solve this problem, then we get an 
approximation in the form (JT3J); cr° is given by (|14|) and cr 1 is obtained from Q17]). 
Since fJ34j) is valid, it follows that the right-hand side of (fTTjl has the form 



m-1 

C 

oc 



M {z=Zj} 



Scr rz rdr. 
r 



Owing to the fact that (JT7j) is linear, the perturbation cr 1 is a linear combination of 
m — 1 solutions of type 



5 Solution techniques 

5.1 Approximate analytical solution 

In this section we construct an analytical solution for the first perturbation term 
cr 1 . Consider a pair of stress functions ((p,ib), such that condition (j3J) is satisfied 

9 (rep) , 3 2 ip 13^ 

a r =cp, a e = — h-^-y-, cr z = — — , a TZ = - — . (35) 

or oz z r or r oz 

These functions wer e constru c ted fr om compatibility equations (jHJ), (|TU|l using the 
standard technique (Washizu, 1982j ). This choice of stress functions simplifies the 



application of variational methods. 



We use the Kantorovich method (see, for example, iKantorovich and Krvlov 



( 1958h ) to reduce the 2-D variational problem to 1-D variational problem. Suppose 



that for ((p,ih), that define the solution cr 1 of (JEj), the following static hypothesis 
is valid 

cp(r,z) = cp 1 (r)cp 2 (z), ^(r,z) = ^iW^fz), (36) 
where <Pi(t),iJ)i(t) are given Kantorovich trial functions, such that 

= epi (r ) = ^i (n) = i\>i (r D ) = 0. 

We use the Kantorovich method to obtain a system of differential equations and 
boundary conditions for (pzi 2 -), tyii 2 -)- This method gives a projection of cr 1 on the 
subspace, defined by (J31j|) . It is more convenient to solve (JT7J) in the form - 
(|32J). For the sake of brevity we describe only the application of the Kantorovich 
method for equation (ffijj) . For equation this procedure can be arranged in a 
similar manner. 
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(37) 



Since cpi(r) and ibi(r) are fixed, the variation of (|36|) gives 
Sep = cpi(r)6cp 2 (z), 6\|> =ih(r)6M> 2 (z). 
Let us rewrite ([27)1 in terms of cp2,ip2 

a 

u T — (h.)-u z Jrdr, 

r . r dz r dr 

^(0)=0, ^(0)=0. 
dz dz 

Here cf( ^2,1^2) and a(5cp 2 , 6i|) 2 ) are defined by (J20J), ((HE), and (l3T |. Using Gauss 
theorem and the fundamental lemma of the calculus of variation ( Washizu . 1982j ) 
we obtain compatibility equations and boundary conditions 



<Pi(r)Ci(e(<p 2 ,^)2))dT = 0, 



\Mr)C 2 (e((p 2 ,ilJ2))dr = 0, 



(3f 



(2s rz -^^)^ 1 (r)dr = 0, 
oz 



(2£rz -^l-^^ lW dr = ) (39) 



3z 



dr 



{z=0} 



{z=H} 



(£e-— )ip 1 (r)dr = 0. 
r 



(40) 



{z=h} 



Here e = C cr and Ci(e), Cz[t) are defined by (fTUjl . (J5QJ). After integration, system 
(J38|) has the form 

d 2 iL> 2 d 2 cp2 d 2 i|>2 d 4 \[>2 

a 1 (p 2 + a 2 ^2 + a 3 -^ r = 0, b 1 (p 2 + b 2 -g 2 - +b 3 \b2 + b 4 -^ 2 - + b 5-^- = 0. (41) 



Elimination of cp 2 from (J4*T|) gives 

k 1 ^ 2 + k 2 



d 2 ib 2 , d 4 Tl> 2 



dz 2 



dz 4 



(42) 



After integration in (jBUjl . pOjl. we obtain 



dij>2 . d 2 ib 2 . d>2 un , n 
dz dz z dz^ 



d^ 2 . d 2 Ti> 2 . d ^2w, x 



du. 



dz 



dz 2 



dz 3 



dr 



■^i(r)dT, 



d 2 ^ 2 



(9i^2 + 92-^- )(H) 



{z=h} 



— tJ>i (r)dr. 



{z=H} 
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Dealing with (|28jh (J2U)) in a similar fashion we obtain compatibility equations and 
boundary conditions for ^>i{z) in (h, H). Finally, taking into account (|5U|). (pTT|). we 
have eight boundary conditions 

^2(0) =0, ^(H)=0, 
dz dz 

dz 

(ei— — + e 2 -r^ + e 3 — -^(0) = 0, 
dz dz z dz- 3 

, d\p2 d 2 \|^2 d 3 i()2,„,, 
(d-p + e 2 —Y + e 3 -r^ H = 0- 
dz dz z dz- 3 

dz dz z dz- 3 



r i , d ^i 



-yihMdr. 



Solution t|>2 (z) is uniquely defined by these equations combined with two ordinary 
differential equations (ODE) (|42|) of fourth order (one equation in (0, fi) and another 
in(h,H)). 

5.2 Numerical solution technique 

In this section we solve problem (JT7j) numerically with the help of the Ritz method. 
Let Gi, i = 1 ,2, N be a system of stresses, satisfying (jlj) and (fl~2|) . Suppose 

N 

CT 1 = Y_ C i°"i- 

i=1 

Then the unknown constants Ct are defined from a system of N linear algebraic 
equations 

dL° 

-r~ lo-=<x°< o" 1 >< o"i >= -h^cr ) < ffi >, i = 1.....N. 
da 

We use the system ffi, i = 1,2,..., produced by means of (J3SJ)- The complete system 
of (p,4> is given by combinations of trigonometric functions 

i - _l n ) cos (^j-j-j-J , i= 1,2,...,N r , j = 0,1,...,N Z , (43) 
t - _ r l 7T ) co H2H/ ' 1 = 1 ' 2 '-' Nr ' ( 44 ) 



1— — ^-ttJ sin ^— J , i= 1,2,...,N r , 



(45) 
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ip = sin fn- ynj cos (k— J , n = 1,2, ...,N r , k = 0, 1 , ...,N Z . (46) 

In order to approximate the stress-jump (JHHj) at the interface {z = h} it could be 
useful to consider additionally discontinuous functions 

<p = sin ( i r - - 7i j H(z), i = 1,2,...,N T , (47) 



\p = sin ( u— — —n | < 



2h 



, z < h 



+ 



zH 



+ 



-HH 



2(H-h) H-h 2(H-h) 



z > h 



,n = 1,2,...,N T . 



(4? 



Here H(z) is the Heaviside function. 

If we put N r = 1 , N z — > oo, then the solution tends to the approximate analytical 
solution from the previous section. 



6 Results of calculations 

The specimen dimensions, applied loads, and material constants are given by H = 8, 
h. = 0.5, T| = 1 , r D = 2, p = 1 , n = 3. 

6.1 Solution of the linear auxiliary problem 

First we compare the approximate analytical solution of (fTTj) (section 5.1) with the 
numerical solution which was obtained using the Ritz method (section 5.2). 

We used (pi (r) = 4h (r) = sin( T ^ r ^ 7t) in to construct the analytical solution. 
Values of constants at, bi, kt, ei, gt are given in Appendix. The numerical solution 
by the Ritz method is performed witn N r = N z = 25 in (J4*3j) — (|48|). 

We use r = 1 , 5 to illustrate on fig. |U the distribution of stress components 
along the z axis. Along r = 1 , 5 the cr z stress is negligibly small for both solutions, 
therefore we do not plot this component. 

Hypothesis (J36|) imposes essential restrictions on the class of solutions, neverthe- 
less the value of shear stress o" TZ is captured very well. The error for the hoop stress 
<Jq and the radial cr T stress is significant in the vicinity of the interface {z = h}. 

The approximate analytical solution gives only a smooth part of the stress field. 
This solution can be more useful in case of smooth changing material properties. 

6.2 Solution of the nonlinear problem 

We investigate the error of the perturbation method o(s) = cr(s) — [cx° + so -1 ). 

Numerical solutions cr(s) of the system (JTJ) — (jHJ) are obtained with the help of 
ANSYS finite element code for a set of s. The geometry of the half of the pipe was 
represented by 3200 axisymmetric PLANE183 finite elements (fig. EJ). We used a 
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A 




Figure 4: Comparison of the approximate analytical solution (dashed line) with 
numerical solution, obtained by the Ritz method (solid line) 




pressure p 



Figure 5: FE mesh and boundary conditions 
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uniform mesh with 160 elements along the axis direction z and 20 elements along the 
radial direction r. This type of element can model creep behavior but the Norton's 
constitutive law is available only as a secondary creep equation. That is, the total 
strain is a sum of elastic strain and creep strain 



If the applied load remains the same with time t, then we have 

e ps e cr as t — > oo, 

and we obtain the steady-creep solution as t — > oo. In calculations we used E = 100 
for Young modulus and v = 0.3 for Poisson ratio to simulate the elastic material. We 
consider the solution to be close enough to the asymptotic solution at the moment 
T = 1 00 of time t. 

The solution cr 1 of the auxiliary problem is given by the Ritz method (see section 
5.2). Here we used series (gSD — (BED with N r = N z = 25. 

The leading term in the asymptotic series is a good approximation even if the 
"small" parameter s equals 0.5 (fig. EJ). With subsequent increase of s the error 
grows dramatically. 



7 Discussion 

The application of the perturbation method to the steady-state creep problem was 
investigated. High performance of this method in prediction of creep response was 
validated. The perturbation method allow one to reduce an initial nonlinear problem 
to the sequence of simpler ones. 

This technique is especially attractive if the unperturbed solution is given in 
closed form as it was in this paper. Another example is the creep response of the 



thick- walled homogeneous pipe under plane stress conditions ([Malininl . 11981( 1 . Such 
solution could be used for perturbation analysis of creep in open-ended pipes. 

The error of the perturbation method becomes substantial when the creep prop- 
erties differ from one another by one order of magnitude. Nevertheless we note that 
asymptotic expansion (|15[) gives a good simplified model of structure response. This 
model treats changes in parameter distribution as jumps of displacements in linear 
elastic material. 

In that way, the solution for every complicated case of parameter distribution is 
represented combination of simple solutions. 
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Figure 6: Inconsistency between o -0 + scr 1 (dashed line) and ct(s) (solid line) 
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Appendix 



Assume cpi (r) = ih (r) = sin( T r n), X\ = 1, r Q = 2. After integration in (|58j) . 
(|40|) we have the following values of constants 

^ w 41 .134, a 2 ^ 3.770, a 3 « -0.237, 

bi « 3.770, b 2 « -0.237, b 3 « 3.948, b 4 « 0.780, b 5 « 1 .778, 
ki « 3.602, k 2 « -0.736, k 3 « 1 .777, 
e 2 = 0, e 3 « -1,777, g 2 « 1,777. 
The ODE (|42jl has four linearly independent solutions 

= e Re ( A ^sin(Im(A)z), ^ = e Re(A)z cos(Im(A)z), 

ty™ = e - Re(A)z sin(Im(A)z), ib 2 v = e- Re(A)z cos(Im(A)z), 
where A is one of solutions of characteristic equation 

k 3 A 4 + k 2 A 2 + kt = 0. 

To be definite, we use 

A = 0.903 + 0.780 t. 
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